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Abstract. A dimer of bound atoms cannot melt, only dissociate. Bulk metals show a well defined first order 
transition between their solid and liquid phases. The appearance of the melting transition is explored for 
increasing clusters sizes via the signatures in the specific heat and the root mean square of the bond 
lengths Sb (Berry parameter) by means of Monte-Carlo simulations of Al clusters modelled by Gupta 
potentials. Clear signatures of a melting transition appear for A*' ~ 6 atoms. Closed-shell effects are 
shown for clusters with up to 56 atoms. The melting transition is compared in detail with the dissociation 
transition, which induces a second and possibly much larger local maximum in the specific heat at higher 
temperatures. Larger clusters are shown to fragment into dimers and trimers, which in turn dissociate at 
higher temperatures. 



PACS. 61.46.-I-W Nanoscale materials: clusters, nanoparticles. 
Thermal properties of small particles, nanocrystals, nanotubes 
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1 Introduction 



The properties of small metal clusters have enjoyed a large 
interest over the past years. Their technical application 
in catalysis stems from the large surface-to-volume ratio 
while their properties differ from those of the bulk mate- 
rial raising the fundamental question about the statistical 
mechanics of finite systems. 

The melting process of small clusters has early on been 
identified as the onset of isomer fluctuations More 
recent investigations on Nii3_a; alloy clusters in the mi- 
crocanonical ensemble ^ show the relation between iso- 
mer fluctuations and the increase in entropy across the 
melting transition. The onset of the melting transition 
is marked by the fluctuations into the lowest energy iso- 
mer configurations which are measured by the root mean 
square bond length fluctuations W sometimes referred to 
as the Berry parameter 5 . The phase space occupied by 
those fluctuations is small for small clusters resulting in 
a maximum in the specific heat at somewhat higher tem- 
peratures. The maximum in the specific heat in turn is 
determined by the onset of isomer fluctuations occupy- 
ing a sufficiently large phase space fraction. A detailed 
overview of the increase of phase space with increasing 
particle number and the classification of isomers in terms 
of potential energy surfaces is given in Ref. 6| . Molecular 
Dynamics (MD) investigations [3 of the melting transi- 
tion of larger Auat clusters with 100 < N < 1000 show 
that the bulk limit is gradually attained in agreement with 
experimental |H] findings. 



Recently the empirical investigation of the melting of 
small Sn pi and Ga ^U] clusters has revealed a possible 
stability of the solid phase of the particles beyond the 
melting temperature of the bulk material. This result was 
interpreted as a consequence of the specific rigid ground 
state structure of the clusters and found support in micro- 
canonical MD calculations for C, Si, Ge, and Sn clusters 
[TTj as well as for isokinetic MD investigations of Snio 
particles [T^ . 

In metals the contribution from the conduction elec- 
trons to the binding energy has to be modelled by many- 
body potentials 13 ,14 , which are numerically more in- 
volved than the thoroughly investigated Lennard- Jones 
systems [51ll5lll6| . A prominent example is the Gupta po- 
tential (GP) [13], which can be derived in the second mo- 
ment approximation from a tight binding model 17 and 
which correctly describes the surface contraction observed 
in metals: 
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Here N is the number of atoms, i and j are atom labels, 
Tij — Tij/rQ — 1, and rij — {fi — rj\ is the modulus of 
the distance between two atoms at positions and r^. 
The parameters have been determined by fitting the ex- 
perimental bulk lattice parameters and elastic moduli ^Sl 
as A = 0.1221 eV, ^ = 1.316 eV, p = 8.612, and q = 2.516 
for Al. Distances are measured in units of the bulk first 
neighbour distance t-q — 2.864 A. 
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The present paper aims to shed light on how the melt- 
ing transition evolves in the limit of small clusters. The 
method is a MC simulation in the canonical ensemble. A 
standard Metropolis algorithm is employed jlHII2()j with 
an update after each random displacement of an atom 
within an interval [0, dmax] in all spatial dimensions, dmax 
is set to yield an MC acceptance rate of 50 to 60 %. The 
resulting temperature dependence is roughly dmax cx 
The boundary conditions arc imposed by a hard wall cube 
with linear dimension L. Runs are performed with sam- 
phng rates (SR) of up to 8 x 10'' steps per temperature 
and atom. The fluctuations on the curves shown in the 
paper are a measure of the statistical error and appear 
near phase transition because of the usual critical slow- 
ing down. The ground state energies and configurations 
obtained within this method are in good agreement with 
earlier results pil21| . The ground state configurations of 
the clusters discussed herein have the same symmetries as 
those of the 9-6 Sutton-Chen potentials \2'2\ . 

An observable commonly studied in the context of 
melting transitions is the Berry parameter ^ 
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where the brackets denote the thermodynamic average in 
the canonical ensemble. The parameter Eq. |j2Jl measures 
the root mean square of the distance between two atoms 
averaged over all pairs. Even short isomer fluctuations 
with a subsequent return to the ground state can leave 
the cluster reordered, i.e., a previously nearest neighbour 
pair Vij may become a second- or third-neighbour pair af- 
ter the fluctuation. Such a reordering leads to a notable 
increase in Sb allowing for the clear identification of a 
melting transition. Note that the Lindemann criterion of 
melting, which measures the atomic fluctuations with re- 
spect to their equilibrium positions, is usually employed 
for bulk systems but is less well suited for cluster systems 

M 

The second observable of interest is the specific heat, 
which can be obtained from the thermodynamic averages 
of the potential energy V and its square: 
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Since the atoms are treated as classical particles the con- 
tribution from the kinetic energy is Ckin = 3/2fcB per 
atom. 

In order to obtain a more intuitive understanding of 
the melting process a real-time visualisation of the sim- 
ulation has been implemented. Figure ^ illustrates snap- 
shots of an Ali3 cluster in a volume of (6ro)'^ at temper- 
atures (a) T = 0, (b) kBT = 0.1 eV, and (c) fceT = 0.4 
eV corresponding to a solid, liquid, and dissociated clus- 
ter, respectively. The graphs show the corresponding nor- 
malised pair distribution functions g{r) (arbitrary units). 
The ground state is icosahedral which is different but close 
to the slightly distorted icosahedral configuration obtained 
by first principle calculations |23|- Contacts of the clusters 




Fig. 1. Normalised pair distribution functions and snapshots 
of an Ahs cluster in a volume of (6ro)^ at temperatures (a) 
kBT = 10"" eV, (b) ksT = 0.09 eV, and (c) k^T = 0.27 
eV corresponding to a solid, liquid, and dissociated cluster, 
respectively. The pair distribution functions in panel (b) and 
(c) are enhanced by a factor of 20 for better visibility. 



with the walls are rare events and the pressure is negli- 
gible in the solid and liquid phases for sufficiently large 
volumes, e.g., L > 4ro for N — 13. 



2 Appearance of the melting transition 

A dimer of bound atoms cannot melt, only dissociate. Bulk 
metals show a well defined first order transition between 
their solid and liquid phases. This phase transition is ac- 
companied by a divergence in the temperature dependence 
of the specific heat indicating the increase in entropy and 
the associated latent heat as well as by a discontinuous 
jump in the Berry parameter Eq. |5J). Figure El shows how 
both signatures evolve [specific heat (a) and Berry param- 
eter (b)] as the particle number is increased from 2 to 10 
atoms. 

AI5 is the smallest cluster with inequivalent bonds in 
its ground state configuration resulting in the abrupt in- 
crease in once isomer fiuctuations set in at around 
ksT ~ 0.013 eV or T ~ 150 K. For A\n with N >6 the 
increase in entropy is sufficiently large around a specific 
temperature to lead to a local maximum in the specific 
heat. As observed previously ^ , the isomer fiuctuations 
leading to the jump in Sb occur at lower temperatures 
than the maximum in the specific heat. The temperature 
of the discontinuity in Sb depends on the energy barrier 
between the ground state configuration and the lowest en- 
ergy isomers [H]. 

Note that the presence of well defined signatures of 
the melting transition in the specific heat for clusters with 
iV < 10 is somewhat unexpected since the experimentally 
investigated Najv clusters appear to show no feature 
in the caloric curves for N < 100 and for Lennard- Jones 
clusters jK] the corresponding signature disappears for 
N < 30. 

It is remarkable that the absence of signatures of a well 
defined melting transition does not imply that the clus- 
ters remain solid up to higher temperatures. For example, 
the real time visualisation of an AI4 particle (within the 
MC/GP approach, see Fig.^for an illustration) at T = 0.1 
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and the Berry parameter (b) of AIn clusters with N = Fig. 3. Temperature dependence of the melting transition for 
2, 3, 4, 5, 6, 7, 10. A clear signature of a melting transition is clusters of size N — 12, 13, 14 in (a)+(c) and N = 54, 55, 56 in 
observed for N > 6 in the specific heat and for N > 5 in the (b)+(d). (a) and (b): specific heat, (c) and (d): Berry param- 
Sb- SR: 4 X 10^ steps per temperature and atom. eter. Volume /vg = 6^. SR: 4 to 8 x 10^ stops [(a)+(c)] and 

1 to 2 X lO'^ steps [(b)+(d)] per temperature and atom. 



eV//cB ~ 1160 K reveals that the cluster fluctuates out of 
its tetrahedral ground state into almost planar configura- 
tions and certainly cannot be considered solid. The small 
size of the phase space is responsible for the absence of 
a signature of the melting in the specific heat while the 
equivalence of all bonds in the ground state configuration 
result in a featureless 6b- 

In this sense all clusters investigated here with > 4 
are in a liquid state already at temperatures below the 
bulk melting temperatures of Al with Tbuik = 933 K ~ 
0.08 eV/fce- For = 2 it is not possible to distinguish 
between a liquid and a solid phase. For = 3 the trian- 
gular structure is stable against fluctuations into coUinear 
isomers up to ksT w 0.15 eV. 

The critical value of the Berry parameter at the melt- 
ing transition |S] was determined as 6b ^ 0.03 — 0.05, 
which is referred to as the modified Lindemann criterion. 
These numbers are consistent with the present observation 
for clusters with N > 5, where at the melting transition 
a jump occurs from values of 6b ^ 0.03 — 0.05 to values 
of 6b > 1.5 JL2,,3,20,25_. For smaller clusters with = 3 
and 4 the temperature dependence of the Berry parameter 
is featureless. For AT = 4 at T = 0.1 eV/fce and for A^ = 3 
at T = 0.15 eV/fce the real time rendering of these clus- 
ters shows that they cannot be considered solid any more. 
At those temperatures the Berry parameter is 6b ~ 0.1. 
These findings suggest that a more general sufficient crite- 
rion for clusters of all sizes not to be considered solid any 
more is 6b > 0.1. The latter is close to the value given for 
the Lindemann criterion 

Multi-step melting |25ll26j and isomer fluctuations 
involving reordered atomic arrangements in the cluster are 
also consistent with that criterion since in both cases at 



least a group of atoms does not remain located at their 
ground state positions when 6b > 0.1. 

Both specific heat and Berry parameter in Fig. ^ 
show an increase for temperatures above ~ 1800 K ac- 
companied by an increase in fiuctuations due to statistical 
errors. As will be discussed further below in detail these 
are signatures of the dissociation transition. 

3 Closed-shell effects 

The infiuence of closed shells on the cohesive energies of 
metal clusters |22ll23ll21|l?i) and their melting points [HI 
I25II20) has been a focus of research for quite some time. 
A closed-shell cluster has a large gap to the first excited 
isomer while adding or removing an atom leads to a num- 
ber of degenerate ground state configurations separated 
by a potential barrier. This manifests itself in a smoother 
specific heat anomaly as well as in a jump in the Berry pa- 
rameter at much lower temperatures as compared to the 
closed-shell counterpart. This is shown in Fig. Eta) and 

(c) for the sets of A^ = 12, 13, 14 and in Fig.|2|;b) and (d) 
for the sets of A^ = 54, 55, 56 atoms. The upper panels (a) 
and (b) show the specific heat, the lower panels (c) and 

(d) the Berry parameter. 

Notably the 14 and 56 atom clusters have a very low 
barrier between isomers. The real time visualisation re- 
veals that these fiuctuations occur in AI14 not only by 
jumps of the 14**^ atom on the surface of the AI13 icosa- 
hedron but by absorption of the 14**^ atom into the outer 
shell and simultaneous pushing of another atom onto the 
surface. In Alsg the 56*'' atom is absorbed into the outer 
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shell even in the ground state configuration (see Sutton- 
Chen 9-6 in Ref. ^|). This leads to the large jump in 6b at 
very low temperatures [dash-dotted lines in Fig. |3fc) and 
(d)] . In Ali4 the phase space is sufficiently large to induce 
an anomaly in the specific heat [dash-dotted line in Fig. 
inja)] at the same temperature. Together with the main 
maximum at higher temperatures this may be referred to 
as a two-step melting mechanism ^25-26-27; . 

The narrower specific heat anomaly at the melting 
transition and the smaller discrepancy between the maxi- 
mum of the specific heat and the jump in the Berry param- 
eter for the larger systems in Fig. OJ^b) and (d) illustrate 
how the thermodynamic limit is gradually approached as 
the cluster size is increased jTJEJIEHl- Note that for larger 
clusters less sampling steps per temperature and atom are 
required to obtain smooth curves. The larger phase space 
of the larger clusters [HI results in the better convergence 
of the observables. 

The canonical ensemble as shown for AI13 in Fig.|2Ia) 
and (b) yields a somewhat lower specific heat |5S| and 
an onset of the isomer fluctuations at lower temperatures 
as compared with the results for the microcanonical en- 
semble in Ref. • The discrepancies can be attributed to 
the energy fluctuations in the canonical ensemble that al- 
low the potential barriers between different isomers to be 
overcome at lower temperatures. 



4 Ideal gas limit 

The well defined high temperature limit, where the system 
has the properties of an ideal gas, yields a test of the 
numerical methods and is essential for determining the 
latent heat of the evaporation transition in Sec. El Figure 
01 illustrates the large T behaviour of the specific heat 
(full line, right scale) and the Berry parameter (broken 
line, left scale) for AI2 (a), AI7 (b), and AI13 (c). Since the 
Boltzmann weight can be expanded as exp{— ^/(fcBT)} = 
1 — V/{kBT) + 0(r~^), the specific heat attains the limit 
as C = I fee + a + 0(r~^) while the Berry parameter 
is = Soo + bT-^ + 0{T-^) with 6^ = 0.3768(1). Here a 
and b are volume and particle number specific constants. 
An analytical scaling analysis reveals that Soa depends on 
the container geometry but neither on its volume nor the 
particle number. 



5 Dissociation 

Between the low temperature liquid phase and the high 
temperature ideal gas limit Fig. 01 shows clear maxima 
both in the specific heat and the Berry parameter for AI2 , 
AI7, as well as for AI13. The feature is generic for all cluster 
sizes and can be associated with the dissociation transi- 
tion. The dissociation anomaly in the specific heat stems 
from the increase in entropy across the dissociation tran- 
sition while the Berry parameter is enhanced through the 
short-time elongation and return of an atom from and to 
the cluster. The latter involves an energy fluctuation and 
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Fig. 4. Temperature dependence of the specific heat (full line, 
left scale) and the Berry parameter (broken line, right scale) 
of a 2 (a), 7 (b), and 13 (c) atom Al cluster. The dissociation 
anomaly is always present while the signature of the melting 
transition evolves with system size. SR: 4 x 10^ [(a) + (b)] and 
2 X 10^ (c) steps per temperature and atom. 



is consequently suppressed in a microcanonical or isoki- 
netic ensemble. 

Figure jsl shows the container size dependence of the 
dissociation transition of an AI13 cluster for volumes of 
L3 = (4ro)3, (6ro)3, (lOro)^, (15ro)3, and (20ro)3. For 

— (4ro)^ the density N/L'^ — 0.0144 mole/cm^ is only 
a factor 7 smaller than that of bulk Al with 0.1 mole/cm'^. 
This leads essentially to a suppression of the evapora- 
tion transition, which is replaced by a smooth crossover 
(see also Fig. Note that even macroscopic particles 
do not exhibit a sharp evaporation transition in a finite, 
constant volume imposing a finite, temperature depen- 
dent gas-phase partial pressure. As a consequence there 
is a liquid- vapour coexistence region which is given by the 
width of the specific heat anomalies shown in Fig. [3 The 
snapshot shown in Fig.EJb) in SecOis taken in the liquid- 
vapour coexistence region of AI13 for = (20ro)'^, where 
liquid fragments of the cluster coexist with evaporated 
single atoms. 

Both the anomaly in the specific heat (upper panel in 
Fig.O and the Berry parameter (lower panel in Fig.jSJ in- 
crease with increasing volume and become narrower. In an 
infinite volume the partial pressure of the gas phase is zero 
and both quantities are expected to exhibit a sharp peak 
at the transition. Note that the fluctuations of the graphs 
in Fig. 13 increase with increasing volume as a consequence 
of the enlarged phase space and thus limit the simulations 
to small volumes. The melting transition is independent 
of the container volume. This is expected since the pres- 
sure in the system is essentially zero when all atoms are 
condensed jSH- 
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Fig. 5. Sample volume dependence of the dissociation transi- 
tion for four different sampling volumes of /r^ = 4^, 6^, 10^, 
15'^, 20^. Both the anomaly in the specific heat (upper panel) 
and the Berry parameter (lower panel) increase with increas- 
ing volume as a consequence of the increased phase space. SR: 
2 X 10'' (L'^/r^ = 4^ 6^ W^) and 4 x lO' (L^/ro = 15^ 20=*j 
steps per temperature and atom. 



For sufficiently small densities or, cquivalently, large 
volumes the latent heat of the melting transition is much 
smaller than the energy released in the evaporation transi- 
tion |2ni- For Ali3 this becomes apparent from the caloric 
curves of the total enclosed system shown in Fig.|H| which 
are obtained by simply integrating the specific heat (Fig.jSJ 
for a fixed volume, i.e., Ev{T) = C{T')dT' - Eb- The 
binding energy per atom for = 13 is Eb = 2.60088 eV. 
Since the transitions are smeared out the determination of 
the latent heat is somewhat ambiguous. Extrapolation of 
the linear segments of the caloric curves above and below 
the melting transition (dashed lines in Fig. ^ and reading 
off the energy difference at the temperature of the specific 
heat maximum yields AE^^it ~ 0.08 eV. 

The evaporation transition is very broad for the vol- 
umes (densities) under investigation and the ideal gas 
hmit (EidcaiiT) = 2/3fcBT, dash-dotted hue in Fig. EJ 
is only attained for fceT > 1 eV even for L/ro = 20. 
For L/tq = 4 the ideal gas limit is reached only for tem- 
peratures much larger than the binding energy per atom, 
i.e., /cbT ::§> Eb- Consequently the total latent heat of the 
evaporation can only be given approximately as Z\_Eovap ~ 
2.26 eV, which corresponds satisfactorily to the potential 
energy expectation value (y)(Tovap) — 2.35 eV near the 
onset of the evaporation transition at /ceTovap = 0.15 eV. 

Results for AI55 (not shown) are similar with slightly 
narrower evaporation anomalies compared with those of 
AI13 with comparable densities N/L^. Closed shell effects 
do not play an observable role in the evaporation transi- 
tion. 
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Fig. 6. Caloric curves as obtained from integrating the specific 
heat data in Fig.\^for sampling volumes of /tq = A^, 6^, 10'^, 
15^, 20'^. The inset shows the determination of the latent heat 
AEn^cit ~ 0.08 eV of the melting transition. The dash-dotted 
line is the ideal gas case reached only for temperatures much 
larger than the binding energy, i.e., k^T S> 2.60088 eV. 



From the above discussion follows that the dissociation 
anomalies in the specific heat usually dominate in size over 
those of the melting transition especially for clusters that 
do not have closed-shell structures. An example is AI7 as 
shown in the middle panel of Fig.0] The jump in 5b is yet 
well pronounced confirming its sensitivity to the melting 
transition [S]. 



6 Fragmentation 

On closer examination the curves of the specific heat of 
AI13 for sufhciently large volumes show another broad lo- 
cal maximum at higher temperatures. For L/r^ — 20 as 
shown in Fig.|SJa) this maximum is found near UbT ^ 0.36 
eV. This effect is readily explained by the fragmenta- 
tion of the cluster into dimers and trimers at tempera- 
tures above the onset of the evaporation transition near 
^sTovap = 0.15 eV. The dimers and trimers in turn disso- 
ciate at temperatures near the second local maximum. For 
smaller volumes (larger densities) the effect is not visible 
because the signatures become very broad. 

To quantify the effect, Fig.^a) shows the maxima of 
the evaporation anomaly of the specific heat as a function 
of the cluster size at constant densities of p = N/L^ = 
1.625 X IO^'^Tq^ , which correspond for A^ = 13 to L/tq — 
20 as shown in Fig.|SJa). The temperature of the specific 
heat maximum increases with increasing size for N > 6 
but for A^ < 6 shows a non-monotonous decrease with in- 
creasing cluster size. In other words, especially the dimers 
and trimers are more stable with respect to larger clus- 
ters. The superposition of their specific heat maxima at 
A^Bl^nax = 0.415(5) eV and kBT^e.. = 0.275(5) eV for 
the dimer and trimer, respectively, leads to the second 
maximum of the curve for L/ro — 20 in Fig. Ela). Fig. 
Ufb) shows a snapshot of AI13 at /cbT = 0.25 eV with 
L/tq — 20, where a dissociated single atom, two dimers 
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Fig. 7. (a) Temperature of the maximum of the specific heat 
at the evaporation anomaly as a function of the cluster size at 
a constant density of p — 1.625 x 10~^r|^^. The line is a guide 
to the eye. The minimum at N — 6 implies that those clusters 
dissociate first while trimers and dimers are stable to much 
higher temperatures. Larger clusters dissociate predominantly 
into dimers and trimers at intermediate temperatures as shown 
by the snapshot in panel (b) for Aha at ksT = 0.25 eV. SR 
for the determination of Tmax-' 10* steps per temperature and 
atom. 

and a trimers are visible together with a (liquid) pen- 
tamere. Qualitatively comparable results are found for 
AI55 (not shown) with feBTmax = 0.23(1) eV at similar 
densities. 

Note that the stabihty of the dimers and trimers to- 
wards dissociation is an entropic effect since the binding 
energy per atom decreases monotonously with decreasing 
size for < 13 [^Uj. Consequently the observed presence 
of dimers and trimers is not in contradiction with the 
results from density functional theory calculations j3ij . 
which find the single atom emission to be energetically the 
dominant dissociation channel. The dimers and trimers 
are at least in part formed by fusion of single evaporated 
atoms. On the other hand, this fragmentation behaviour is 
not directly transferable to experiments on Al clusters be- 
cause the Gupta potentials do not yield the correct planar 
structure for A^ < 5 [2SIEI1- 



7 Conclusions 

The results from the Monte-Carlo simulation of AIat clus- 
ters modelled with many-body GPs reveals the appear- 
ance of a distinct feature of the melting transition in the 
specific heat and the Berry parameter for A^ > 6. The 
energy fluctuations in the canonical ensemble lead to an 
onset of isomer fluctuations at lower temperatures than 
in the microcanonical ensemble. AIn clusters with A^ > 4 
are liquid at the bulk melting temperature. The present 
analysis suggests that a generalised sufficient criterion for 
clusters not to be considered solid any more is a Berry 
parameter of i5b > 0.1. Larger clusters with closed-shell 
configurations exhibit sharper signatures of the melting 
transition than others. 

For higher temperatures clusters of all sizes undergo a 
dissociation transition which is accompanied by container 
volume dependent anomalies both in the specific heat and 
Berry parameter. Dimers and trimers are more stable to- 
wards dissociation than larger clusters. The details of the 



features depend on the potential but the qualitative re- 
sults are generic. For example, MC simulations with GPs 
and parameters for Au yield qualitatively very similar re- 
sults. 

The author thanks M. Blom, P. Schmitteckert, G. Schnei- 
der, D. Schooss, M. Vojta, and P. Wolfle for instructive 
discussions. The work was supported by the Center for 
Functional Nanostructures of the Deutsche Forschungsge- 
meinschaft within project D1.5. 
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